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Abstract 

The uncertainty of the theoretical prediction of the B — > X s ^f branching ratio at 
NLL level is dominated by the charm mass renormalization scheme ambiguity. In 
this paper we calculate those NNLL terms which are related to the renormalization 
of m c , in order to get an estimate of the corresponding uncertainty at the NNLL 
level. We find that these terms significantly reduce (by typically a factor of two) the 
error on BR(Z? — > X s ^f) induced by the definition of m c . Taking into account the 
experimental accuracy of around 10% and the future prospects of the B factories, we 
conclude that a NNLL calculation would increase the sensitivity of the observable 
B — > X s ^f to possible new degrees of freedom beyond the SM significantly. 
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1 Introduction 



The branching ratio of B — > X s ^ is a very sensitive probe for new degrees of freedom be- 
yond the standard model (SM) (for a review, see [1]). Within supersymmetric extensions 
of the SM for example, one can derive stringent bounds on the parameter space of these 
models [2-8]. Clearly, such bounds will be most valuable when the general nature of the 
new physics beyond the SM will be identified at the forthcoming LHC experiments. 

Because of the heavy mass expansion that is valid for inclusive decay modes, the 
decay rate of B — > X^ is dominated by the perturbatively calculable partonic decay rate 
r(6 — > X s j). QCD corrections to the latter, due to hard-gluon exchange, are the most 
important perturbative contributions; they were calculated in the past up to the next-to- 
leading logarithmic (NLL) level [9-18]. Subsequently, also electroweak corrections were 
calculated [19-22]. After completion of these computations, it was generally believed that 
the theoretical uncertainty of the branching ratio is below 10%. 

However, as first pointed out in 2001 in [23], there is an additional uncertainty in the 
NLL results for T(b — > X s j) which is related to the definition (renormalization scheme) 
of the charm quark mass. Technically, the charm quark mass depencence enters through 
the matrix elements (sjlOi^b) which in the context of a NLL have to be calculated up 
to 0(a s ). As these matrix elements vanish at the lowest order, the charm quark m c only 
enters (through the ratio m c /m&) at 0(a s ). As a consequence, the charm quark mass 
does not get renormalized in a NLL calculation, which means that the symbol m c can 
be identified with m CjPO i e or with the MS mass m c (/i c ) at some scale /i c or with some 
other definition of m c . Formally, all these assignements are equivalent, as they lead to 
differences which are of order a 2 s . 

Note that in contrast to the c-quark mass the 6-quark mass does get renormalized in 
a NLL calculation and we choose to express all the following results in terms of m^ po i e . 
In this respect we do not follow ref. [23], where the m^is mass was used. Unless stated 
otherwise, the symbol rrib stands for mb iPO i e in all the formulas in this paper. Numerically, 
we use rrib = 4.8 GeV throughout. 

Numerically, it turns out that the NLL result for F(b — > X^) strongly depends on 
which mass definition of the charm quark mass is used in the NLL expressions. To 
illustrate this, we first identify m c with m CjPO i e as it was done in all analyses before the 
paper of Gambino and Misiak [23]. Numerically, we use m c ^ \ c /mb tP o\ e = 0.29 which is 
based on the mass difference m^ po i e — m CiPO i e = 3.4 GeV fixed through the heavy mass 
expansion of and nip and mb po i e = 4.8 GeV. The corresponding branching ratio then 
reads [23] 

BR[5 -> X sl } Ei>mb/20 = 3.35 xlO" 4 . (1) 

As the charm quarks which are propagating in a loop have a typical virtuality of mb/2, 
the authors of Ref. [23] suggested to use m c ( / u c ) with /i c e [m c ,m b ] instead of m C]PO i e . A 
typical value for the corresponding ratio is m c (/i c ) / m fepo i e = 0.22. Using this value, the 
branching ratio gets increased w.r.t. JT} by about 11% [23]: 

BR[B ^ X sl ] Ei>mb/20 = 3.73 xl0~ 4 . (2) 

In a recent theoretical update of the NLL prediction of this branching ratio, the uncer- 
tainty related to the definition of m c was taken into account by varying m c jnib in the 
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conservative range 0.18 < m c /mb < 0.31 which covers both, the pole mass (with its 
numerical error) value and the running mass m c (/x c ) value with \i c e [m ci m^\ [24]: 

BR[B -> X s7 ] = (3.70 ± 0.35| mc/m6 ± 0.02| C km ± 0.25| param . ± 0.15| scalc ) x 10~ 4 . (3) 

There exists a large number of measurements of the inclusive decay B — > X s 7 [25-30] 
and the present experimental accuracy has reached the 10% level [31]: 

BR[B -> X s7 ] = (3.52 ± 0.30) x 10" 4 . (4) 

In the near future, more precise data on this mode are expected from the B factories. 
Thus, it is mandatory to reduce the present theoretical uncertainty accordingly. A sys- 
tematic improvement certainly consists in performing a complete NNLL calculation . This 
is, however, a very complicated task (for discussion and some results see [32-35]) and a 
certain motivation is needed to enter such an enterprise. In the present paper we try to 
give such a motivation: By calculating those NNLL terms which are induced by renormal- 
izing the charm quark mass in the NLL expressions, i.e. those terms which are sensitive 
to the definition of the charm quark mass, we show that the large error at the NLL level 
related to the m c definition gets significantly reduced. As this error is the dominant one 
at the NLL level (see eq. (j3J)), we conclude that a complete NNLL calculation will drasti- 
cally improve the theoretical prediction of the branching ratio. We stress here that in the 
present paper we only make a statement about the reduction of the error at the NNLL 
level, and not about the central value of the branching ratio; this remains the topic of a 
complete NNLL calculation! 

The remainder of this paper is organized as follows. In section |21 we discuss in some 
detail how to calculate the NNLL terms induced by renormalizing m c in the NLL results. 
In order to make the paper self-contained, we first list in section |3] the structure of the 
NNL results and then we present the analytical results for the new terms discussed in 
section |21 Finally, in section 01 we numerically investigate by how much the error related 
to the definition of m c gets reduced at the NNLL level. 



2 NNLL terms related to m c renormalization 

As already explained in the introduction, the matrix elements M"2 t (m c ) = (s7|0i i 2(/x&)|6) 
only start at order 0(a 4 ), or, in other words at the NLL order 2 . As a consequence, the 
definition of m c is not fixed at this order, because m c does not get renormalized. This is 
also true for the bremsstrahlung contributions Mj )r 2 ems (m c ) = (sjg\Oi t 2(lJ>b)\b) , which are 
needed up to 0(g s ) for a NLL calculation. In this section we concentrate on the virtual 
terms M^lrric), as the extension to the bremsstrahlung contributions M^2 ms { m c) is 
straightforward. 

2 In the present paper we use the operator basis as first introduced in ref. [11]. /i& denotes the renor- 
malization scale of O(mft). 
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Figure 1: Left frame: Typical 5m c insertion diagram (see text). Right frame: Typical 
diagram with a charm quark self-energy insertion (see text). 



When going to NNLL precision, the matrix elements Ml l ^{m^) are needed to 0{a 2 s ). 
At this level, there are - among many other diagrams - counterterm contributions to 
these matrix elements, induced by the renormalization of m c (see the left frame of Fig. 
[TJ. The complete set of such diagrams is generated by inserting the operator —i5m c ip c ip c 
in the 0(a s ) diagrams of in all possible ways. The sum (JMj'f'^mc) ■ 5m c of all these 

insertions can be obtained by replacing m c — > m c + 5m c in the 0(al) results 2 {m c ) , 
followed by expanding in Sm c up to linear order: 

M$ t(£) K + Sm c ) = M™ t(t \m c ) + 8M^ {e \m c ) ■ 5m c + 0((5m c ) 2 ) . (5) 

As 5m c is ultraviolet divergent, the matrix elements (m c ) are needed in our appli- 

cation up to order e 1 , as indicated by the notation in eq. (J5J). 

The explicit shift 5m c depends of course on the renormalization scheme. When aiming 
at expressing the results for (mc) i n terms of rh c (fib), the shift reads (Cf = 4/3) 

5m c (fi b ) = - "^^ C F -m c (fi b ) . 

47r e 

On the other hand, when the result is expressed in terms of m CiPO i e , the shift reads 



om CjPO i e = C F - + 3 In — + 4 m c , po i e 

Air \ e mr 1 



c 



The infinities induced by the 1/e terms in Sm c get cancelled in a full NNLL calculation, 
in particular by self-energy diagrams depicted in the right frame of Fig. ^ As we do not 
perform a full NNLL calculation, we suggest to consider self-energy insertions, where the 
self-energy S(p 2 ) is replaced by £i(p 2 = m 2 ). 

The Ex-part of the self-energy Ex(p 2 ) is defined through the decomposition of the full 
unrenormalized self-energy S(p 2 ) as 

S(p 2 ) = S 2 (p 2 )(^-m c ) + S 1 (p 2 ). (6) 
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At the one-loop level, the corresponding pieces £f and of the renormalized self-energy 



~*R 



arc 



^(p 2 ) = £ 2 (p 2 ) + 5Z C , Sf (p 2 ) = ^(p 2 ) + 5m, 



R(J2\ 



(7) 



where Z c = 1 + 5Z C denotes the wave function renormalization constant of the charm 
quark. Eq. (JJJ) implies that the sigularities in ^[^^(mc) • 5m c cancel when combined 
with the diagrams with Si(j» 2 ) insertions. However, for general p 2 , the function Hi(p 2 ) 
depends on the gauge parameter £: 



4n 



-CpfTlc I 3 



e mi 



?7l 2 \ / p 2 

p 2 I m 2 



As for p 2 = m 2 the self-energy piece £i(p 2 = m 2 ) is gauge- independent, we add £i(p 2 = 
m 2 ) insertions to 8Mj^^ (m c ) ■ 5m c . 

These momentum independent Si(jo 2 = m 2 ) insertions can be straightforwardly aborbed 
into 5m ( i s insertions: 



5m eS 
°'' t c,pole 

5mf^ b ) 



^(p 2 = m 2 c ) + 5m CtPO i c = , 



E x (p = m c ) + 5m c (/j, b ) 



i4 



47T 



C F 31n^ + 4 m c (// 6 ). 



Finally, if we wish to express the matrix elements M^" (m c ) in terms of m c (/z c ), the 
shift reads 



,ir, . _ a a (ji b ) {., | /' 



(Jmf(Aic) 



47T 



m 



(9) 



3 Analytical results 



Before turning to the contributions induced through the renormalization of the charm 
quark mass, which are NNLL terms, we first summarize the structure of the NLL result 
for the branching ratio for b — > X s ^f. We write the decay width for b — > X s ^ using a 
photon energy cut Eq — ^(1 — $) = E max {l — 5) as 



r(6 -> X s ^) E > Eo = r(6 -> S7) + r(6 -> sjg) E ^> Eo 



(10) 



where the two parts are defined as follows: 



G 



r(6 -> s7^)s 7 >i?o = i^j\ V ts V tb\ 2 ^e m rril vole A 



32tt 



D = C? )eff (/i 6 ) + 



47T 



f^ 1)cfr (^) + E^ 0)eff 

V i=i 



(M6 



^ + 7i7 



(O)eff ln 



^6, 



16 



4 



A = ( e -*WM^MWM _ i) |cf eff (// 6 )| 2 + £ ^MG^M&W • 

(12) 



i 



The expressions for the Wilson coefficients Cj(/^,) can be found in [36]. Their numerical 
values we take from table 5.1 in ref. [37]. Writing the results in this specific form, the 
functions fij(S) and are understood to be taken from [11] and not from the original 
paper [10] where the results were parametrized differently. 

Following common practice, we write the branching ratio (without taking into account 
non-perturbative corrections) as 

BR(b -> x sl ) Ei> _ Eo = r( r ( 7^x!5r BRT ' (13) 

where the semileptonic decay rate is given by 

= (14) 



1927T 3 _ 

g(z) = 1 — 8 z + 8 z 3 — z 4 — 12 z 2 \n(z) is the phase-space factor and the function 

2a s (m b ) f(z) 



with 



3tt g(z) 



2 x /25 239 25 2 \ , . , / 4 , 17 , 

f(z) = -(1 - z 2 ) - — z + — z 2 j + z \n(z) ( [ 20 + 90z--z 2 + -z 3 

+z 2 ln 2 (z) (36 + z 2 ) + (1 - z 2 ) ( ^ - ^ z + ^ z 2 ) ln(l - z) 

z)-tt 2 



3 3 3 

-4 (1 + 30 z 2 + z 4 ) \n(z) ln(l - z) - (1 + 16 z 2 + z 4 ) (6 Li( 



-32z 3/2 (l + 



7T 



2 



4 Li( v^) + 4 v/i) - 2 ln(z) In 



,1 + x/i, 



accounts for 0(a s ) QCD corrections. We note that m c is understood to be the pole mass 
in eq. IjUjl . 

We now turn to that part of NNLL corrections which is responsible for the reduction 
of the charm quark mass renormalization scheme dependence, as explained in section EJ 
We first turn to terms SMj 1 ^^ induced by m c renormalization in the matrix elements 

Mi 1 .^ 6 . To this end, we need up to oder e 1 . In [10] have calculated these matrix 

elements up to terms e°, using Mellin-Barnes representations for generalized propogator 
to obtain analytic results in the form of the series in z = (m c /mb) and L = ln(z). As in 
these calculations the expansion in e was the last step, it is straightforward to calculate 
U P f° order e 1 . 

In order to get finite results for these matrix elements, we add counterterms related to 
operator mixing as in ref. [10], adapted however, to the operator basis defined in ref. [11]. 
This step leads to ,ren , which we decompose as in ref. [10]: 

= (sj\O r \b) £ 1„ 15 - Z£ eh > 15 - U In ^rf + rf + r« £ ) . (15) 
47r V 81 fib 81 fib fMb ' ) 



We obtain for r 2 = r^ + er 2 (note that r\ = — | r-z): 

rf = - — (-48 - 3 L 2 - L 3 + 5 vr 2 + 9 L (-4 + vr 2 ) + 36 C(3)) z + — vr 2 z 3 / 2 

2 243 27 V v ; sv V 27 

8 4 
+ — (L 3 - 6 L(-2 + 7t 2 ) + 18 + 2 vr 2 ) z 2 - — (9 - 182 L + 126 L 2 + 14 vr 2 ) z 3 



i 2 7T 



' U 1 2 (-15 - 3 L - 3 L 2 + vr 2 ) z + 2 (-3 L 2 + vr 2 ) z 2 + - (-7 + 3 L) z 3 



27 L 3 V ) V ) 3 

(16) 

r?> = + ^ vr 2 - A (-18180 + 75 L 4 + 3240 vr 2 + 46 tt 4 - 30 L 2 (-24 + 7 vr 2 ) 

+9000 C(3) + 120 L(-66 + 14 vr 2 + 27 C(3))) z - — tt 2 (-49 + 6 L + 24 ln(2)) 2 3/2 

' 81 

2 

H (48 L 3 - 15 L 4 + 24 L 2 (-3 + vr 2 ) - 24 L(3 + 5 ?r 2 ) + 1116 + 36 vr 2 

81 ^ 

+40 vr 4 + 432 C(3)) z 2 - tt 2 z 5/2 + — (22705 - 2484 L 2 + 4536 L 3 - 6036 vr 2 
sv 'I 81 729 V 

+6 L(-1783 + 192 vr 2 ) + 8208 C(3)) z 3 



i Z 7T 



27 



221 



9 



+ (l5 L 2 - 6 L 3 - 4 L(-9 + tt 2 ) + 186 - 10 tt 2 - 36 C(3)) 



—2 (—3 — 6 L + 3 L + 2 7r + 2 Lit 2 + 18 C(3) j z 2 
4 

+ - (-67 + 66 L + 9 L 2 + 12 vr 2 )r 3 
9 



(17) 



In these formulas we retained all terms up to order z 3 , as higher order terms contribute 
much less than 1%. Nevertheless, in the numerical evaluations in section 0] all terms up 
to z 6 were included. 

At the level of the decay width, the implementation of the contribution coming from 
renormalization of the c-quark mass in the virtual contributions is (according to eq. (jHJ)) 
most easily achieved by replacing r 12 in eq. (|llj) by rf\ + ^ r i,2, where 

Ar li2 = 5m c ^(rW + er? 2 ) . (18) 

At the NLL order, the bremsstrahlung corrections to the decay width are encoded in 
the quantities fij(5) (see eq. ffT2^) ). which correspond to the interference terms (Oi,Oj). 
In the following, we calculate the shifts Afy to these quantities induced by the renor- 
malization of the charm quark mass. In principle, we calculate the decay width using a 
photon energy cut 5 = 0.9 (see eq. (JT0J) ) . However, as all bremsstrahlung contributions 
which contain charm quark loops are finite for 8 — > 1, we can approximate these terms 
by putting 5 = 1. Numerically the relative error is of order 10~ 4 . 

We first calculate the shift A/27- To this end, we shift the charm quark mass in 
the matrix element of (s r yg\02\b) as in eq. © and then work out the interference with 
(s^g\0'j\b) . Because of the 1/e term in 5m c , the result is ultraviolet singular. In a 
full NNLL calculation this singularity gets cancelled when combined with self-energy 



6 



insertions in the charm quark lines in the matrix element of 2 - We therefore do the 
phase space integrations involved in the derivation of f 2 i ( or A/27) i n d=4 dimensions. 
As only the matrix element of 2 depends on m c , the shift A/27 can be constructed by first 
considering the quantity /27 itself. Using the integral representation for the building block 
for photon and gluon emission from the c-quark loop [10], one obtains after integration 
over all but one of the phase space parameters 



/ 



27 



I l^L \ I dx dy duu 2 (\ — u)(\ — x)y 
9 \rrib 



2c 



[x(l -x)]- e r(2 + e) e e7+iOT 



uy - 



x(l— x) 



+ fc) 



(19) 



Here x, y are Feynman parameters and u is the remaining phase space parameter, < 
x,y,u < 1. To solve the integrals, we use the Mellin-Barnes representation for the 
generalized propagator 



uy 



x{\ 



X 



IT] 



-1-e 



2inT(l + e) 



ds- 



T{-s)T(l + e + s) 
(uy) l +- 



uyx[l 



x 



appearing in eq. p9|l . 7 denotes the integration path parallel to imaginary axes which 
hits the real axes somewhere between (— 1 — e) and 0. Closing the integration path in the 
right s-half plane, one gets an expansion for f 2 ? in z = [m c jmb) 2 and L = ln(z). 

The shift A/27 is then obained as 



A/ 



27 



5m r 



dz df: 



27 



dm r dz 



d(0) 
27 



+ e /2 d 7 (1) ) 



(20) 



To summarize, the NNLL contributions due to renormalization of m c in the (O2, O7) 
interference are taken into account by replacing / 2 7 — > /^ + A/27 i n eq. (fl^jl. Explicitly, 
we find: 



f (o) 
J27 



+ -(7 - 2rr 2 + 6L + 2L 2 )z + (rr 2 -2L- L 2 )z 2 



12 8 



6tt 2 - 4L + QL 2 )z s 



-(-6 + lOL)^ 4 



+ ^(13 + 70L)z 5 + ^(32 + 63L)z 6 



(21) 



,d(0) 

J27 



+ 



-(13 - 2rr 2 + 10L + 2L 2 ) + 2(-l + vr 2 - 3L - L 2 )z 
9 L8 

-(-37 - 18tt 2 + l%L 2 )z 2 + -(-7 + 20L)z 3 
4 3 



(22) 



+ ^(27 + 70L)^ 4 + i(85 + 126L) 2 5 



fd(l) 

J27 
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— (165 - 52tt 2 - 4(-15 + rr 2 )L - 18L 2 - 8L 3 - 72((3)) 
48 



+ -(2(-12 + rr 2 )L - 3L 2 + 4L 3 + 6(-5 + 4tt 2 + Q((3)))z 



- ^(-7+15vr 2 + (15 + 27r 2 )L-15L 2 + 4L 3 + 36C(3))^ 2 (23) 

- -4(235 + 48tt 2 - 204L + 36L 2 )z 3 
27 v ' 

+ ^(- 10076 - 4200tt 2 + 19425L - 3150L 2 )^ 4 

+ —(-3554 - 4200vr 2 + 20685L - 3150L 2 )/ 
250 v ; 

Note that / 27 in eq. (}2~Tj) is an expanded version in z of the integral expression for / 2 7 i n 
ref. [11]. We further note that / 28 — — I/27, fn = — I/27, As = 7^/27; the same relations 
also hold for the respective A/y (see for instance, [38]). 

Finally, we turn to the shift A/22 related to the (0 2 , (^-interference. To derive 
this quantity, one has to perform the shift m c — > m c + 5m c only in one of the two 
interfering one- loop amplitudes M 2 brems = (s^g\0 2 \b). To this end, one writes integral 
representations for both, M2 rems and dM2 rems / dm c . A/22 is then represented as a five 
dimensional integral (4 Feynman parameters and one phase space parameter), which can 
be solved by double Mellin-Barnes techniques (see for instance [39]). Omitting the detail 
of this calculation, the terms induced by renormalizing m c in the (0 2 , 2 ) bremsstrahlung 
terms are implemented in eq. (|T2"j) by replacing / 22 — > / 2 2 + A/22, where 

A/ 22 = 2^ (24) 
rrib nib v ' 

Explicity, we get: 

f$ = 0.04938272 + (16.64197 + 1.887290L - 0.4444444L 2 - 0.09876543L 3 )^ 
+ (57.92026 + 47.67037L + 1.185185L 2 + 3.134737L 3 + 0.05925926L 5 )^ 2 
+ (-93.12628 + 32.36078L - 12.95977L 2 + 1.777778L 3 - 0.2962963L 4 )/ 
+ (11.92082 - 11.21491L + 2.074074L 2 - 0.5925926L 3 );? 4 (25) 
+ (0.6482797 - 4.160089L + 0.1810700L 2 - 0.329218LL 3 )z 5 
+ (-1.125313 - 4.320604L - 0.2444444L 2 - 0.3456790L 3 )/ , 



/ 22 (0) = 18.52926 + 0.9984013L - 0.7407407L 2 - 0.09876543L 3 

+ (163.5109 + 97.71112L + 11.77458L 2 + 6.269473L 3 + 0.2962963L 4 + 0.1185185L 5 )z 

+ (-247.0180 + 71.16280L - 33.54596L 2 + 4.148148L 3 - 0.8888889L 4 )^ 2 

+ (36.46839 - 40.71149L + 6.518519L 2 - 2.370370L 3 )z 3 (26) 

+ (-0.9186906 - 20.43831L - 0.08230453L 2 - 1.646091L 3 )^ 4 

+ (-11.07248 - 26.41251L - 2.503704L 2 - 2.074074L 3 )^ 5 , 



/ 22 {1) = 41.24600 - 7.794263L - 0.7525535L 2 + 0.3950617L 3 + 0.07407407L 4 
+ (234.4505 + 44.95451L - 64.68047L 2 - 0.5200208L 3 
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- 4.498135L 4 - 0.05925926L 5 - 0.08559671L 6 )z 

+ (-368.8104 + 245.2526L - 95.81857L 2 + 28.84099L 3 - 3.851852L 4 + 0.5728395L 5 )^ 2 

+ (-3.708986 - 93. 16023L + 34.50854L 2 - 7.670782L 3 + 1.283951L 4 )z 3 (27) 

+ (-45.10910 - 54.73258L + 13.39517L 2 - 4.073160L 3 + 0.8779150L 4 U 



4 

+ (-92.60340 - 74.33726L + 10.32328L 2 - 4.786008L 3 + 1.094650L 4 )z 5 . 

We decided to give the expansion coefficients in these equations in numerical form, because 
the exact results are somewhat lenghty. We note that /° 2 in eq. (J23j) is an expanded version 
in z of the integral expression for fi2 in ref. [11]. We further note that /n = and 
fi2 = — 5/22! the same relations also hold for the respective Afij. 

These analytical results are defined parts of the complete NNLL contribution which can 
be used within a future NNLL calculation. 



4 Numerical results 



In the following analysis we show that the NNLL terms, induced through the renormal- 
ization of m c , drastically reduce the error related to the definition of the charm quark 
mass in BR(6 — > X s j). To illustrate this feature as clearly as possible, we take the fixed 
values shown in Table d for the input parameters. In particular, we use the fixed ratio 



m b = 4 


8 GeV 


m CtPO i e /m b 


= 0.29 m z = 91.18 


7 GeV 
















a s {m z ) = 0.119 


"em 


1/137.036 


\v; s v tb /v cb \ 2 = 0.95 


BR S \ = 


= 10.49% 



Table 1: Input parameters used in the numerical analysis. 

m c, P oie/ m b,po\e = 0.29. Furthermore, we always leave the semileptonic decay width, which 
enters the branching ratio for b — > through eq. (fT3*j) . expressed in terms of m cpo i e as 
given in eq. (fH^ . In this way the m c definition dependence of the BR(6 — > X s ^) only 
comes from the numerator in eq. (fT3*|) . For our studies, we neglect electroweak corrections 
and non-perturbative effects. As already mentioned, in the bremsstrahlung contribution 
we use 5 = 0.9 for the lower cut in the photon energy (see eq. (jlOjO . 

Starting from m CjPO i e = 0.29 • 4.8 GeV = 1.392 GeV, we first calculate 7Ti c (m CjPO i e ), 
using the one-loop expression 



&s (^c,pole 



7T 



(28) 



To get m c (/i c ) for an arbitrary scale (typically between 1.25 GeV and 5 GeV), we use 
two-loop running (with 5 flavours) according to 



m c (fj, c ) = m c {no) 



V «»(W))/ 



1 + 



2/3 Q 2PI ) 



4tt 



(29) 



with /i = m C)PO ie- Numerically, we get the values shown in tabled In Figure[2]our results 
are given for three different values of fi b , where [i b represents the usual renormalization 
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m c (1.25)/m 6 = 0.257 



m c (2.5)/m b = 0.214 



m c (5.0)/m 6 = 0.187 



Table 2: m c (/i c )/mb for /z c = 1.25,2.5,5 GeV using m c ^ \ e /mb = 0.29 as input. 

scale of the effective field theory. We compare the branching ratio for b —>■ X s ^ within 
the pole and the MS scheme for the charm quark mass. Within each vertical string the 
solid dot represents the branching ratio using m CjPO i e , while the open symbols correspond 
to rh c (fi c ) for fi c = 1.25 GeV (triangle), /i c = 2.5 GeV (quadrangle) and /i c = 5.0 GeV 
(pentagon), respectively. 
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Figure 2: BR(6 — > X s ~f) for three values of /if,. For each value of fib the left string shows 
the NLL results for m CtPO i e (solid dot) and for m c (/i c ) with \x c = 1.25; 2.5; 5.0 GeV (open 
symbols). The right strings show the corresponding NLL results supplemented by the 



Sm c mass insertions and the £i(p 



insertions (see text for more details). 



For each ^ the left string shows the value of the branching ratio at the NLL level, 
while the right string shows the corresponding value where in addition Sm c mass insertions 



and £i(p 



insertions were taken into account, as explained in detail in section |21 



Because the combination of these insertions is zero by construction for the pole scheme 
(see eq. (jHJ)), the solid dots are at the same place in the left and the right string for a 
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given value of 

From Figure 121 we see that the error related to the charm quark mass definition gets 
significantly reduced when taking into account NNLL terms connected with mass inser- 
tions. Taking as an example the results for fib = 5 GeV, we find that at the NLL level the 
branching ratio evaluated for m c (2.5 GeV) is 12.6% higher than the one based on m cpo i e , 
in agreement with ref. [23]. Including the new contributions, these 12.6% get reduced to 
5.1%. 

A remark concerning the remaining NNLL terms is in order: As these terms give 
contributions to the branching ratio which (up to terms of order a^) do not depend on 
charm quark mass definition, the error related to m c in the full NNLL result is expected 
to stay essentially the same as estimated in the present paper. 

However, to obtain a NNLL prediction for the central value of the branching ratio, it 
is of course necessary to calculate all NNLL terms. 

Summing up, we have shown that the relatively large error related to the definition 
of the charm quark mass in the NLL result for BR(6 — > X s ~/) gets significantly reduced 
(typically by a factor of 2) at the NNLL level. Taking into account the present experi- 
mental accuracy of around 10% and the future prospects of the B factories and also of 
possible Super- B factories [40,41], we conclude that a future NNLL QCD calculation of 
the b — > X s 7 branching ratio will significantly increase the sensitivity of this observable 
to possible new physics. 
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